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ABSTRACT 

I While supernova remnants have been identified as the most likely sources of the 

' galactic cosmic rays, no conclusive observational evidence for this association ex- 

^SJ I ists. We show here that IceCube has the possibility of producing incontrovertible 

evidence by detecting neutrinos produced by the cosmic ray beam interacting 
with the hydrogen in the vicinity of the supernova shock expanding into the 
interstellar medium. We show that the observational information on gamma ray 
fluxes from the Cygnus region, although limited, is sufficient to pinpoint the 
expected event rate of the neutrinos associated with a single source of 0.5 Crab 
at the TeV level to within a factor of two, between 2 and 3.8 neutrinos per year. 
' Finally, we note that recent gamma-ray observations reveal the presence of at 

least three and possibly up to eight such sources, raising the possibility of seeing 
O ' more than 10 neutrinos per year from these sources alone. 

c3 I 1. Cosmic Neutrinos Associated with Galactic Cosmic Rays 

ps| . Cosmic accelerators produce particles with energies in excess of 10*^ TeV; we do not 

^ ! know where or how. The flux of cosmic rays observed at Earth follows a broken power 

^ [ law. The two power laws are separated by a feature dubbed the "knee". Circumstantial 

^ I evidence exists that cosmic rays, up to perhaps EeV energy, originate in galactic 

^ supernova remnants. Any association with our Galaxy disappears in the vicinity 

Q ■ of a second feature in the spectrum referred to as the "ankle". Above the ankle, the 

^ ■ gyroradius of a proton in the galactic magnetic fleld exceeds the size of the Galaxy and 

it is generally assumed that we are witnessing the onset of an extragalactic component 
in the spectrum that extends to energies beyond 100 EeV. While the enigmatic nature 
of the highest energy cosmic rays has been widely advertised, it is also a fact that the 
5^ \ origin of the galactic cosmic rays has not been established. We will show here that the 

positive identiflcation of supernova remnants as the cosmic ray accelerators is possible 
with the kilometer-scale neutrino observatories now in the planning or construction 
phase. 

Simple energetics point at the supernova origin of the galactic cosmic rays. By 
integrating the observed flux we can obtain the energy density Pe of cosmic rays in 
the galaxy from the relation that flux = velocity x density, or 

The answer is that pe ~ 10~^^ ergcm"^. This is also the value of the corresponding 
energy density B'^/Stc of the microgauss magnetic fleld in the galaxy. The acceler- 
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ator power needed to maintain this energy density is 10~^^erg/cm^s given that the 
average containment time of the cosmic rays in our galaxy is 3 x 10^ years. For a 
nominal volume of the galactic disk of 10^^ cm^ this requires an accelerator delivering 
lO^^erg/s. This happens to be 10% of the power produced by supernovae releasing 
Iq5i gj,g gygj-y 3Q years. This coincidence is the basis for the idea that shocks produced 
by supernovae expanding into the interstellar medium are the origin of the galactic 
cosmic rays. 

Can we observe neutrinos pointing back at the accelerators of the galactic cosmic 
rays? The conversion of the 10^° erg of energy into particle acceleration is believed 
to occur by diffusive shock acceleration in the young (1000-10,000 years) remnant 
expanding into the interstellar medium. If high-energy cosmic rays are indeed asso- 
ciated with the remnant, they will interact with hydrogen atoms in the interstellar 
medium to produce pions that decay into photons and neutrinos. These may provide 
us with indirect evidence for cosmic ray acceleration. The observation of these pionic 
gamma rays has been one of the motivations for neutrino as well as ground-based 
TeV-energy astronomy. 

The HESS telescope opened a new era in astronomy by producing the first re- 
solved images of sources in TeV gamma rays, particularly, in this context, of the 
supernova remnant RX J1713.7- The observed flux suggests that HESS may 

have identified the first site where protons are accelerated to energies typical of the 
main component of the galactic cosmic rays. Although the resolved image of the 
source reveals TeV gamma ray emission from the whole supernova remnant, it shows 
a clear increase of the flux in the directions of known molecular clouds. This is sug- 
gestive of protons, shock accelerated in the supernova remnant, interacting with the 
dense clouds to produce neutral pions that are the source of the observed increase 
of the TeV photon signal. The image shows filaments of high magnetic fields con- 
sistent with the requirements for acceleration to the energies observed. Follow-up 
observations of the source in radio- waves and X-rays failed to identify the population 
of electrons required to generate TeV photons by purely electromagnetic processes; 
for a detailed discussion see Ref. E) 

Fitting the observed spectrum by purely electromagnetic processes is challenging 
because the relative height of the inverse Compton and synchrotron peaks requires 
very low values of the 5-field, inconsistent with those required to accelerate the 
electron beam to energies that can accommodate the observation of 100 TeV photons. 
Nevertheless, an exclusively electromagnetic explanation of the non-thermal spectrum 
is certainly not impossible and even favored by some One can, for instance, 
partition the remnant in regions of high and low magnetic fields that are the respective 
sites of acceleration and inverse Compton scattering. 

A similar extended source of TeV gamma rays tracing the density of molecular 
clouds has been identified near the galactic center. Protons apparently accelerated by 
the remnant HESS J1745-290 diffuse through nearby molecular clouds to produce a 



signal of TeV gamma rays that trace the density of the clouds®. Detecting this source 
in neutrinos will be challenging because it is relatively weak — its TeV luminosity is 
only of order 0.1 Crab. This is^resumably due to the larger distance to the source as 
compared to RX J1713.7-3946^ Moreover, its large angular size on the sky results 
in a larger background of atmospheric neutrinos. So far, HESS has not claimed the 
discovery of pionic gamma rays and finding neutrinos as a smoking gun for cosmic 
ray acceleration in supernova remnants remains of interest. 

Supernovae associated with molecular clouds are a common feature of associations 
of OB stars that exist throughout the galactic plane. Although not visible by HESS, 
intriguing evidence has been accumulating for the production of cosmic rays in the 
Cygnus region of the Galactic plane from a variety of experiments: 

• The observation of the Cygnus region by the HEGRA air Cherenkov telescope 
resulted in the serendipitous discovery of a TeV 7-ray soured TeV J2032+4130, 
with an average flux of ~3% of the Crab Nebula^ and a hard spectrum. Es- 
pecially intriguing is its possible association with Cygnus 0B2, a cluster of at 
least 2700 (identified) young, hot stars with a total mass of ~10^ solar masse P. 
At a relatively small distance of approximately 1 kpc, this is the largest massive 
Galactic stellar association. 

• The Whipple Observatorjff^ confirmed an excess in the same direction as 
J2032+4130, although with considerably larger average flux above a peak en- 
ergy response of 0.6 TeV. Their latest results"'""'"' reveal a TeV hot spot (integral 
flux ~8% of the Crab) that is displaced about 9 arc minutes to the northeast of 
the TeV J2032+4130 position. 

• A re-analysi^^ of the radio surveys of the region revealed a weak non-thermal 
shell-like supernova remnant-type object with location and morphology very 
similar to the HEGRA source. It could be the cosmic ray engine that powers 
the OB association. 

• The Milagro Collaboration reports an excess of events from the Cygnus region 
at the 10.9a level^^. The observed flux within a 3° x 3° window centered at the 
HEGRA source is 70% of the Crab at the median detected energy of 12 TeV. 
Such a flux largely exceeds the one reported by the HEGRA Collaboration, 
implying that there could be a population of unresolved TeV 7-ray sources 
within the Cygnus 0B2 association. In fact, they report a hot spot, christened 
MGRO J2019+37, at right ascension = 304.83°±0.14stat±0.3sys and declination 



= 36.83° ± 0.08stat ± 0.25syi^. A fit to a circular 2-dimensional Gaussian 
yields a width of 0.32 ± 0.12 degrees, which for a distance of 1.7 kpc suggests a 
source radius of about 9 pc. The brightest hotspot in the Milagro map of the 
Cygnus region, it represents a flux of 0.5 Crab above 12.5 TeV. Interestingly, 
the Tibet AS-gamma Collaboration has observed a cosmic ray anisotropy from 
the direction of Cygnus, which is consistent with Milagro's measurement''^^ 



As for the HESS sources, the observations suggest the production of cosmic rays 
as well as a variety of opportunities for neutrino production. The model proposecP^ 
for MGRO J2019+37 is that of a cosmic ray beam which escapes from the OB star 
cluster and interacts with a molecular cloud positioned a few degrees to the southeast. 

The discussion emphasizes the importance of observing the neutrinos from the 
decay of charged pions that accompany the TeV photons if these are indeed the 
decay products of neutral pions produced by cosmic rays interacting in the interstellar 
medium in the vicinity of the supernova remnant. The particle physics is simple: After 
oscillations over cosmic distances the neutrino beam consists of equal fluxes of muon, 
electron and tau-neutrinos and their antiparticles. Since proton-proton collisions yield 
two charged pions for every neutral pion, oscillations cause the flux for each neutrino 
flavor to equal one half of the gamma ray flux. Because the protons transfer on 
average 0.2 of their energy to secondary pions, and the four leptons in the charged 
pion decay chain tt e + Ue + v^) + z/^ take roughly equal energy, neutrinos 

with .05 of the cosmic ray energy are produced. Similarly, photons with 0.1 of the 
proton energy are made from the decay of neutral pions. Accelerators producing 
cosmic rays reaching the "knee" must produce photons with energies up to 100 TeV 
and neutrinos up to half that energy. This requirement is consistent with observations 
of RX J1713.7-3946 and MGRO J2019+37 discussed above. 

With 677 optical sensors in place since February 2000, the existing AMANDA 
detector has been collecting neutrinos at a steady rate of four per day. These at- 
mospheric neutrinos are the byproduct of collisions of cosmic rays with atmospheric 
nitrogen and oxygen nuclei in the northern atmosphere. Note that at the South Pole 
one observes neutrinos that originate in the northern hemisphere, using the earth as 
a filter to select neutrinos from other particles. No photons, or any other particles 
besides neutrinos, can traverse the whole planet to reach the detector. The signals 
from the atmospheric neutrinos do not yield information about astronomy yet, but 
they are calculable and have been used to calibrate the detector. As in conventional 
astronomy, AMANDA looks through the atmosphere for cosmic signals and the data 
are scrutinized for hot spots in the northern sky that may signal cosmic sources. The 
fluxes from the sources discussed in this paper are small and, as has been known for 
some time, their detection requires the construction of neutrino detectors of kilometer 
scale. 

Starting in the Antarctic summer 2004-2005, IceCube deployments have been 
steadily augmenting the AMANDA instrumentation. As of January 2007 IceCube 
consists of 1424 digital optical modules distributed over 22 strings and 54 surface 
cosmic ray detectors. The hardware and software have worked out of the box and 
revealed the first atmospheric neutrinos. The instrumented volume of IceCube already 
exceeds one quarter of a kilometer cubed and will deliver a kilometer-square year of 
integrated data by 2008-09 provided the deployments remain on schedule. We will 
show that this represents the first opportunity to observe the muon neutrinos possibly 



accompanying the gamma rays observed in the Cygnus region. IceCube has sensitivity 
to the Southern HESS sources by observing the showers initiated by electron and tau 
neutrinos. A quantitative estimate is not possible as the fraction of these events for 
which the direction can be reconstructed with degree accuracy has not been published. 

With the initial deployment and operation of strings in the Mediterranean, the 
Antares collaboration has observed its first neutrinos and is establishing a technology 
for constructing a kilometer-scale observatory in the Northern hemisphere. It will 
be ideally positioned to observe the supernova sources in the Southern part of the 
galactic plane including RX J1713. 7-3946, the most promising source so far. 

In this paper we will first review the qualitative arguments why the observed 
TeV gamma ray fluxes from the supernova remnants discussed above are consistent 
with those expected from a generic source building the steady galactic flux observed. 
We will argue that although the estimates are very uncertain, observation of the 
accompanying neutrino flux requires kilometer-scale detectors. Uncertainties in the 
calculation are associated with the propagation of the cosmic rays, with the value of 
the magnetic fields and the age of the remnant. Additionally, in the case of MGRO 
J2019+37, the spectral slope has not been measured. After investigating the wide 
parameter space of models for MGRO J2019+37, we show that the neutrino flux can 
be predicted up to a factor of 2 once we match the fluxes at 12.5 TeV to the Milagro 
data and limit the GeV flux by imposing the constraint that EGRET did not observe a 
GeV counterpart This represents our main result. Using the published sensitivity 
of IceCube to muon neutrinoJ^ we predict 2 ~ 4 events per year in a degree angular 
bin centered on MGRO J2019-I-37 on a background of atmospheric neutrinos of 2.5 
events per year. Confirmation of cosmic ray acceleration should emerge after a few 
years of data taking. 

2. Secondary Pionic Gamma and Neutrino Fluxes from 
Supernova Remnants 

A simple estimate is sufficient to conclude that the observed gamma ray flux 
associated with the supernova remnant RX J1713. 7-3946 and MGRO J2019+37 is 
consistent with the energetics required from typical sources of galactic cosmic rays. 
The emissivity in pionic gamma rays (photons produced per cm^ s) resulting from a 
density of accelerated protons Up interacting with a density of hydrogen atoms n in 
the interstellar medium is 
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The emissivity of photons is simply proportional to the density of cosmic rays 
np{> 1 TeV) (~ 4 X 10"^'* cm~^ for energy in excess of 1 TeV) and the target density 
n of hydrogen atoms. The proportionality factor is determined by particle physics: 
{Et^I Ep) ~ 0.2 is the average energy of the secondary pions relative to the cosmic ray 
protons and Xpp = {napp)~^ is the proton interaction length (cTpp ~ 40 mb) in a den- 
sity n of hydrogen atoms. (We here assumed a generic E^^ spectrum of the protons, 
for different spectral indices the quantity {ET,/Ep)_js generalized to the spectrum- 
weighted moments for pion production by nucleons^; see later). 
The total luminosity in gamma rays is given by 

W ,o 1 

LJ>1 TeV) = — ^ 10^^ photons s"^ (5) 
P 

The density of protons from a supernova converting a total kinetic energy W of 
10^° erg to proton acceleration is approximately given by W/ p, where we will assume 
that the density in the remnant is not very different from the ambient energy density 
p ~ 10~^^ergcm~^ of galactic cosmic rays. This approximation is valid for young 
remnants in their Sedov phas^^. 

We thus predict a rate of TeV photons from a supernova at a distance d of 1 kpc of 

in-ii I photons\ I ^ \ f 'IT' \ I d \ 
d{lnE) ^ '~ And? ~ \ cm2s / VlO^Oerg; Ucm-sj l^lkpcj ' ^ ' 

This prediction is credible because the number of TeV photons predicted coincides 
with observations of the supernova remnant RX J1713. 7-3946 by the HESS array of 
atmospheric Cherenkov telescopeJII'. 

Within the precision of the astrophysics it is adequate to assume that muon neu- 
trinos and muon antineutrinos are produced at the level of one half the gamma flux 
as explained in the introduction. We thus anticipate an event rate of ~1.5 detected 
neutrinos per decade of energy per km^ year, a result readily obtained from the rela- 
tion 

dN events , /neutrinosX /A^\ 

where the last factor represents, as before, the probability that the neutrino is de- 
tected. It is approximately 10~^ for the TeV energy considered here. This estimate 
may be somewhat optimistic because we assumed that the sources extend to 100 TeV 
with an spectrum. The number of neutrinos is uncomfortably small and sug- 
gests a more sophisticated calculation that allows for an estimate of the range of 
fluxes allowed within observational constraints. We perform such an analysis for the 
observation of the Cygnus region by IceCube, whose acceptance is simulated following 
Refs.EniandEH 

Alternative candidates have been suggested for the sources of the galactic cosmic 
rays for, instance microquasars. The above argument suggests that they should have 



left their imprint on the Milagro skymap and they did not. It is very suggestive that 
the Milagro sources are the cosmic ray accelerators. 



3. Remnants in the Sedov Phase: Fluid Effects on the Proton Spectrum 

The excess of events detected by the Milagro Collaboration in the Cygnus region 
of the Milky way may be due to a shell-type supernova remnant (SNR). In such an 
object, a spherical shock sweeps through the local ISM, colliding charged particles 
accelerated by Fermi shock acceleration into the ambient medium. These collisions 
lead to the production of neutral and charged pions that decay to photons and leptons, 
respectively, as discussed in Section 2. 

The injected proton spectrum Q{E) will be modified by fluid processes acting on 
the particles as they traverse the volume of the SNR during acceleration. Normally, 
considering only the shock acceleration gives some probability for particle escape 
during each successive traversal of the shock front. We consider two additional mech- 
anisms that allow the protons to escape and thereby modify the final proton spectrum: 
convection and diffusion. Convection is the escape of particles due to currents within 
the object, while diffusion represents the random walk of particles resulting in a net 
motion down the particle density gradient. Since the mean free path scales with 
energy, diffusion will steepen the high-energy tail of the proton spectrum, with high 
energy particles escaping before those with lower energies. Here we assume Bohm 
diffusion, which gives us the dependence on energy of the escape time. 

We take the equation describing the final proton spectrum Np to be; see e.g. 
Ref.El 

- ^^^^ - 7(^7' 

where t{E) is the characteristic escape time of the protons. We find the solution in 
the steady-state limit as 

N{E)=Q{E)t{E), (9) 

where t{E) = rfluid(^) when rfluid(^) < tsNR and t{E) = tgNR when rfiuid(^) > tsNR,; 
tsNR is the age of the supernova remnant and Tn-aidiE) is the escape time due to fluid 
effects. 

Tfluid(-£') is given by the harmonic mean of the convective and diffusive escape 
times 

1 1 1 

(10) 
(11) 
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Tfluid Tdiflf Ta 

where 

kR 



^shock 

and 




Tdiff = (4524 years) x — x - x (12) 



for k the compression factor of the gas, R and v the radius and velocity of the shock, 
and B the magnetic field in the region. The parametrization of Tdiff is from Ref. |23| 
and assumes Bohm diffusion, hence the scaling. We choose and hold constant 
R = 10 pc and Vghock = 1000 km/s. While above approximations may be inadequate 
for a detailed study of the sources, they are sufficient for deriving the neutrino flux 
from the astronomical observations. The calculation is insensitive to the detailed 
astronomical parameters as we will show further on. 

The result of modifying the initial proton spectrum by these mechanisms is equiv- 
alent to introducing a gradual cutoff at the energy where the behaviour of the 
Bohm diffusion takes hold-essentially, is superimposed on the injected spectrum. 
In the limits we describe above, this transition energy is determined by the intersec- 
tion of the fluid escape time curve with the age curve of the SNR (Figured]). We 
modify the spectrum by multiplying it by whichever of the two times represents a 
smaller escape time at each energy. This means that if Tfiuid > ^snr for a given en- 
ergy, the spectrum is simply multiplied by some constant to be determined by the 
normalization and its shape at that energy is unmodified from the original injection 
spectrum. While convective transport is included in the model for completeness, the 
characteristic time for convective escape does not depend on energy. Therefore, vary- 
ing the convection parameters as well as the age of the remnant is redundant as their 
effects on the final gamma-ray spectrum are almost identical up to normalization. 
Hence the holding constant of the convection parameters. 
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Figure 1: Comparison of SNR age (blue dash-dotted), fluid escape time (solid black), diffusion time 
(red dotted) and convection time (green dashed) for a sample parameter set of age — 1000 years, 
_B— field = 1/iG, Vghock = 1000 km/s and R = lOpc . In our model we modify the spectrum in 
proportion to whichever of the age or the fluid escape time is lower at a particular energy. 



This model introduces a number of parameters whose precise values for the Cygnus 
object are unknown although there is a region of the parameter space considered 
most physically likely: radii of order several parsecs, velocities of order several tens 
of thousands of km/s, magnetic fields of order between several and several tens of 
microGauss, and ages between 500 and 10,000 years. However, since we are normaliz- 
ing the resulting 7-ray spectrum to the observations made in ~10TeV 7-rays by the 
Milagro experiment and constraining the spectrum to be below the fluxes measured 
by the EGRET satellite experiment in the GeV range, the model becomes highly 
degenerate. There are essentially two regimes: 

1. Diffusion begins to affect the spectrum at energies higher than the Milagro 
observation at 12 TeV. This happens when the 5-field is large or the age of the 
object is small. In this case, the loss of particles will happen at an energy at 
which the flux was already very low, resulting in little effect on the neutrino 
event rates from an undiffused spectrum; 

2. Diffusion begins to affect the spectrum at energies lower than the Milagro point. 
This means that the Milagro observation is of a flux already reduced by this fluid 
mechanism, resulting in a potentially large boost to the flux in the GeV range. 
However, the largest boosts to the midrange energies can only be obtained by 
considering potentially unphysical values of the fluid parameters — B of order 
several tenths or hundreths of a microGauss, for example. This range is also 
constrained by the EGRET measurements. 

4. Determining the Gamma-Ray Spectrum from the Proton Spectrum 

In this section we discuss a model of the proton acceleration that allows us to derive 
a pionic gamma-ray flux. We assume a shock accelerated proton input population 
that can be described by a power law with an exponential cutoff. To accommodate 
the "knee" in the cosmic-ray spectrum, we put the cutoff at 1 PeV: 



The choice of cutoff, although chosen to be at the same scale as the knee at 3 PeV, 
is not critical to the following analysis. The flux is so low at PeV energies that there 
are zero events regardless of the precise location of the cutoff. 
The basic formula for the gamma spectrum is 



where Np is the proton spectrum modified by fluid effects (previous section) and F 
is a fragmentation function defined in Ref. |2H describing the energy transfer of the 
accelerated protons to neutral mesons (pions and etas) and subsequently, to 7-rays. 



Qp = AE'e-^'^ 



(13) 




(14) 



We plotted the gamma-ray flux predicted by the above model for a range of 
parameters, normalized to the Milagro observation at 12.5 TeV and constrained by 
the EGRET observations in the ~100 MeV-10 GeV energy range. Following Ref. [5] 
the EGRET source 3EG J202 1+3716 provides the lower-energy constraint. It is clear 
that 3EG J2021+3716 is either the counterpart to MGRO J2019+37, or, if it is 
not, the true counterpart must have a smaller flux than it in order to have escaped 
detection. In either case, the EGRET measurement gives us a clear upper limit in 
the low-energy 7-ray range. 

We took three sample input proton spectra s = —2.0, —2.2, —2.34 and plotted 
each for four values of the magnetic field ranging from 0.1 to 50 fiG. The age of the 
SNR was then chosen to be the value that maximized the gamma-ray flux for the 
previously chosen set of parameters, subject to the EGRET constraint (Figures [2H3]). 
For s > —2.34, there were no gamma-ray spectra that did not exceed the EGRET 
upper limit. Indeed, s = —2.34 (Figure HI) proved such a constraint on the fluid 
parameter space that the variation in the allowed spectra produces only a very small 
variation in neutrino event rates in IceCube. Due to the degeneracy of the model 
discussed earlier, we did not need to vary the variables associated with convection. 




E, (TeV) 

Figure 2: 7-ray spectra with injection s = —2.0. All lines are for 5000 year old SNR. The black 
dashed line is for a magnetic field of 0.1 /iG, the red dotted line is for a magnetic field of 1 fiG, the 
green dash-dotted line is for a magnetic field of 10 /iG, and the blue solid line is for a magnetic field 
of 50 fiG. The inverted triangles are show the EGRET measurement of 3EG J2021-I-3716 and the 
solid box is the Milagro measurement of MGRO J2019-I-37. 
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Figure 3: 7-ray spectra with injection s = —2.2. The black dashed hne is for a magnetic field of 
0.1 and age of 1400 years, the red dotted line is for a magnetic field of 1 fiG and age of 5000 years, 
the green dash-dotted line is for a magnetic field of 10 /iG and age of 5000 years, and the blue solid 
line is for a magnetic field of 50 fiG and age of 5000 years. 
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Figure 4: 7-ray spectra with injection s ~ —2.34. The black dashed line is for a magnetic field of 
0.1 /iG and age of 100 years, the red dotted line is for a magnetic field of 1 fiG and age of 1000 years, 
the green dash-dotted line is for a magnetic field of 10 fiG and age of 1000 years, and the blue solid 
line is for a magnetic field of 50 /iG and age of 1000 years. 



5. Neutrino Event Rates 



We calculated the neutrino flux dN^/dE^, from the gamma flux by the method of 
Ref. [251 assuming the two are equal up to some constant that also includes the effects 
of oscillations (see Ref. 1^ for an alternative method and results). We can neglect the 
energy dependence of the oscillations due to the large distance to the source. Given a 
neutrino flux, the event rate in IceCube is simply the convolution of the flux with the 
energy-dependent area and the probability of an incident neutrino inducing a muon 
with a visible track. The area is taken to be the effective area after quality cuts 
on the IceCube data referred to as level 2 cuts in Ref. EH 



The event rates in IceCube are calculated in this way are shown in Figures [5H71 
The rates are within the range 2 < dN/dt < 3.8 events per year with the IceCube 
threshold at 50 GeV. This is in large part due to the fact that the Milagro observa- 
tion strongly constrains the flux in the energy range 1-20 TeV, where the neutrino 
detection probability is highest, resulting in similar predictions for dissimilar SNR 
characteristics. The irreducible atmospheric background, due to neutrinos produced 
in the Northern atmosphere in cosmic ray showers, is calculated using the results of 
Ref. [27] and is found to be approximately 2.5 neutrinos/year at 50 GeV. Hence in 15 
years of operation we predict 4.9 cr < N/ y/ A'^atmo < 9.3 a and if the higher end of 
the predicted event rate range is realized, 5 a is possible in 4.3 years. 

In closing, we note that the Milagro collaboration^ has recently detected mul- 
tiple additional sources besides MGRO J2019+37, most with approximately equal 
fluxes of 0.5 Crab. The sources with possible counterparts in the GeV range indicate 
a spectral index of ~ — 2.3, bearing out our observation in Section 4 that spectral 
indices steeper than —2.34 exceed the EGRET signal from the possible GeV coun- 
terpart to MGRO J2019-I-37 for all model parameters. If we compute the flux of 
neutrinos from the Milagro sources (not including the Crab Nebula) detected with 
post-trial signiflcance of greater than flve a assuming a power-law index of —2.34, we 
get a total event rate in IceCube of 6.9 neutrinos/year. If we also include the more 
tentative sources, the event rate increases to 11.5 neutrinos/year. In the long run, a 
correlation analysis of the IceCube and Milagro skymaps should make the detection 
of these sources likely. 
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Figure 5: Events due to the gamma-ray spectra with injection a = —2.0 shown in Fig. [21 
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Figure 6: Events due to the gamma-ray spectra with injection a = —2.2 shown in Fig. [S] 
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Figure 7: Events due to the gamma-ray spectra with injection a ~ —2.34 shown in Fig. ID 
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